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Summary. We study the quantum summation (QS) algorithm of Brassard, H0yer, 
Mosca and Tapp, see [1], which approximates the arithmetic mean of a Boolean 
function defined on A*' elements. We present sharp error bounds of the QS algorithm 
in the worst-average setting with the average performance measured in the L, norm, 
q G [l,cxj]. 

We prove that the QS algorithm with Af quantum queries, M < N, has the 
worst-average error bounds of the form ©(In M/M) for q = 1, &{M~^^'') for q € 
(l,oo), and is equal to 1 for q = oo. We also discuss the asymptotic constants of 
these estimates. 

We improve the error bounds by using the QS algorithm with repetitions. Using 
the number of repetitions which is independent of M and linearly dependent on q, 
we get the error bound of order for any q £ [1, co). Since f2{M~^) is a lower 

bound on the worst-average error of any quantum algorithm with M queries, the 
QS algorithm with repetitions is optimal in the worst-average setting. 



1.1 Introduction 

The quantum summation (QS) algorithm of Brassard, H0yer, Mosca and 
Tapp computes an approximation to the arithmetic mean of the values of a 
Boolean function defined on a set of iV = 2" elements. An overview of the QS 
algorithm and its basic properties is presented in the first two sections of [4] . 
In Section 1.2 we remind the reader of the facts concerning the QS algorithm 
that are needed in this paper. 

* The research of the second and third coauthors was supported in part by the Na- 
tional Science Foundation (NSF) and by the Defense Advanced Research Agency 
(DARPA) and Air Force Research Laboratory under agreement F30602-01-2- 
0523. 
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The QS algorithm enjoys many optimaHty properties and has many ap- 
phcations. It is used for the summation of real numbers which in turn is an 
essential part for many continuous problems such as multivariate and path in- 
tegration, and multivariate approximation. The knowledge of the complexity 
of the quantum summation problem allows us to determine the quantum com- 
plexity of many continuous problems, such as those mentioned above, see [6] 
and a recent survey [3]. 

The QS algorithm has been studied in the two error settings so far: 

• worst-probabilistic in [1, 4], 

• average-probabilistic in [4] . 

These settings are defined by taking the worst case/average performance with 
respect to all Boolean functions and the probabilistic performance with respect 
to outcomes of the QS algorithm. 

It turns out that the QS algorithm is optimal in these two settings. The 
corresponding lower bounds for the Boolean summation problem were shown 
in [5] for the worst-probabilistic setting, and in [7] for the average-probabilistic 
setting. In particular, we know that the QS algorithm with M quantum 
queries, M < N, has the error bound of order in the worst-probabilistic 
setting. 

In this paper we study the worst-average setting. In this setting, we 
take the worst case performance over all Boolean functions and the average 
performance over all outcomes of the QS algorithm. The average performance 
is measured in the Lg norm, q G [l,oo]. This setting is analogous to the 
randomized (Monte Carlo) setting used for algorithms on a classical computer. 
The worst-average setting also seems to be quite natural for the analysis of 
quantum algorithms. 

As we shall see, the results depend on the choice of q. Obviously, for 
larger q, the effect of the average behavior becomes less significant. In fact, 
the limiting case, g = oo, leads to the deterministic case (modulo sets of 
measure zero). Not surprisingly, for g = oo, the results are negative. 

In what follows we indicate error bounds for large M. Since we always 
assume that M < N, this means that for M tending to infinity we also let 
N tend to infinity. To make error bounds independent of N, we take the 
suprcmum over A'^ > Af in the corresponding definitions of the errors. When 
we speak about the sharpness of error bounds, we usually take a large M and 
select a still larger N and a Boolean function for which the presented error 
bound is sharp. 

The worst-average error e^°'^~'^^^{M) of the QS algorithm with M quan- 
tum queries satisfies: 




Furthermore, the asymp- 



totic constant is 2/7r for M — 2 divisible by 4. 
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• For g e (1,00), we have e™'"^^s(Af) = —1- ]■ Furthermore, the 

asymptotic constant is approximately ( JJ^ sin'^~^(x)dx/7r) for M — 2 
divisible by 4 and q close to 1. 

• For g = 00, we have e^^-^^^{M) = 1. 

The error bounds of the QS algorithm arc improved by the use of repeti- 
tions. Namely, we repeat the QS algorithm 2n -h 1 times and take the median 
of the outputs obtained as the final output. This procedure boosts the success 
probability of the approximation at the expense of the number of quantum 
queries. We show that with n independent of M and linearly dependent on g, 
we decrease the QS algorithm error to be of order . Hence, the use of 
repetitions is particularly essential for large q since we change the error bound 
0(M~^/*) without repetitions to the error bound 0{M~^) with repetitions. 
The constant in the last big O notation is absolute and does not depend on q 
and M. 

The error bound of order is optimal. This follows from the use of, for 
instance, Chebyshev's inequality and the fact that the lower bound i7(M~^) is 
sharp in the worst-probabilistic setting, see also [7]. Hence, the QS algorithm 
with repetitions is optimal in the worst-average setting. 



1.2 Quantum Summation Algorithm 

The quantum summation QS algorithm of Brassard, H0yer, Mosca and Tapp, 
see [1], approximates the mean 



1 



i=0 

of a Boolean function / : {0, 1, . . . , A''— 1} {0, 1}. Without loss of generality 
we assume that A'' is a power of two. 

The QS algorithm uses M — 1 quantum queries. The only interesting case 
is when M is much smaller than N . The QS algorithm returns an index 
j e {0, 1, . . . , M - 1} with probability 

^Va^{MQas) ( . _2 f - ^af)\ , . _2 / + CTa^ ) 



2M2 V \ M J \ M 

see [4] for the detailed analysis of the QS algorithm. Here 

M 

6a f = arcsin ^/aj and ctoj, = — Oaf 
We will also be using 

Saf = min{ \aaf] - (Ja^aa, - [aafl }. 
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Clearly, s^^ € [0, i] and s^^ = iff a^j, is an integer. We shall usually drop 
the subscript / and denote Oa = , (Ja = <^af i = when / is clear from 
the context. 

Knowing the index j, we compute the output 

a/U) = sm ' 



(1.1) 



on a classical computer. The error is then given by 



sm — — — sm ' 



M \ M 



As in [4], we let l-i{-,f) denote the measure on the set of all possible out- 
comes of the QS algorithm which is defined as 

M(A/) = ^p/(j) VAc{0,l,...,M-l}. 

Let Am denote the set of all possible outputs of the QS algorithm with M—1 
queries, i.e.. 

Am = I sm^ (^] : j = 0, 1, . . . , M - ll. 



M 



Let 



p/(a)=M(|ie{0,l,...,M-l}:sin2(^^^ ""}'^) ^ 

denote the probability of the output a. Note that a = sin'^(7rj/M) 
= sin^(7r(M - j)/M). Hence if j ^ and j ^ M/2 then pf{a) = pf{j) + 
PjiM-j). 

In what follows we let Bjv denote the set of all Boolean functions defined 
on {0,1,..., AT - 1}. 



1.3 Performance Analysis 

The error of the QS algorithm in the worst-probabilistic and average-probab- 
ilistic settings has been analyzed in [1, 4]. In this paper we analyze the error of 
the QS algorithm in the worst-average setting. This corresponds to the worst 
case performance with respect to all Boolean functions from Bjv and the 
average performance with respect to all outcomes. This average performance 
is mcasiued by the expectation in the Lq norm, q G [l.oo], with respect 
to the probability measure of all outcomes provided by the QS algorithm. 
As mentioned before, we make the worst-average error independent of N by 
taking the supremum over N > M. That is, the worst-average error is defined 
as: 
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for q € [1, oo), 

1/9 



M-1 



e— 8(M)= sup max !«/ " a/OOl^ 

N>M ' ^ — ' 



• for g = oo, 



'(M) = sup max max |a/ -a/(j)|. 

N>M /e»Nj:p/O)>0 



It is easy to check that for q = oo, the QS algorithm behaves badly. Indeed, 
if M is odd, we can take / with all values one, and then a/ = l,pf{Q) = 
and a/(0) = 0. Hence e^°'-'*^s(M) = 1. If M is even, we take / with only one 
value equal to 1, and then a/ = ^/N , p/(M/2) > and a/(M/2) = 1. Hence, 
|a/ - a/(M/2)| = 1 - 1/iV and e™'-'^^s(M) = 1. 

That is why in the rest of the paper we consider q G [l,oo). As we shall 
see the cases q > 1 and q = 1 will require a different analysis and lead to quite 
different results. 

1.3.1 Local Average Error 

We analyze the local average error for a fixed function / G]BArforl<(j'<oo, 

el^^{f,M)=(Y,PfU)\af-af{jW] = p/(a)|a/ - a|M . 

(1.2) 

We first analyze the case g' > 1. 

Theorem 1. Let q e (l,oo). Denote a = af. If aa e Z then e^^s(/,M) = 0. 
If aa ^ Z then 

-irs^/M 



e-«(/,M)«-^^^i-P^ / " sin(:r)«-2|sin(x + 2^,)|«da; < 

n_Lon A ^^ 7^'^~^sin(7^Sa) , sin^(7rSa) f C ■q-2(^^^ 

(l + 2(l-(5g,2)) + — — I 2(l-(5,,2) + gy^ sm« {x)dx 



(1.3) 

with S„ = \_aa\ - (Ja and Sa = (Ta - {(Ta] 



* Note that the last integral is finite. It is obvious for q > 2. For q G (1,2), the 
only singularities are at the boundary points and are of the form x''~^ for x 



approaching 0. The function a' ^ is integrable since q > 1. 
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Proof. If (Tj, e Z then it is shown in [4] that there exists a G Am such that 

a = Uf and pf{a) — Sa^af for all a € Am- Then e^^{f, M) = as claimed. 

Assume that CTo ^ Z. Using the form of Pf{j) from Section 1.2, we 
rewrite (1.2) as 



(erW.M))'=E=^ 



j=0 



+ 



sm 



7r(i - tJa) 



M 



sm 



We have 

M-l 



E 

j=0 



sm 



7^0' + 
M 

M 

= E 



9-2 



sm 



7r(j + CTg) 

M 



71" (j - o-q ) 



9-2 



9-2 



sm 



sm 



+ cra) 



M 



M 



sm 



7r(M - j + o-g) 
M 



9-2 



sm 



7r(M - j - gg) 
M 



Using the 7r-periodicity of | sin a; |, we see that the last sum is equal to 



M 

E 



sm 



M 



sm 

M-l 

E 

j=0 



M 



sm 



M 



9-2 



sm 



+ (Jg) 

M 



Therefore 



with 



(1.4) 



M-l 



sin 



- o-g^ 



M 



9-2 



We split ^M.g as 



■ 



M-l 

= E 

j=0 



M 



sin 



sm 



7r(j + o-g) 



M 



9-2 



9-2 



sm 



TTJ 

M 



M 



Observe that ^S^j^^ is the rectangle formula for approximating the integral 
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/ 



sin(a; — 6'a)|'' ^Ism^x + BaWdx. 



aai/M,77laa]/M] 



The error of the rectangle quadrature for fc G N and an absohitely con- 
tinuous function / : [o, 6] ^ M whose first derivative belongs to Li([o, 6]) 
satisfies 



Thus defining h{x) = | sin(a; — 0o)|' ^|sin(x + 0o)|'^ and Da = [0, tt] \ 
[tt [aa\ /M, TT [cTa] /M] and using the error formula above for the subintervals 
[0, TT [aa\ /M) and (tt \aa\ /M, tt], we get 



< 



b — a 



\f'{x)\dx. (1.5) 



— S' 



/ h{x)dx <— \h'{x)\dx. 

J Da ^ J Da 



Define H{x) = h{x + (9„) = | sin(.T)|'' ^| sin(x + 26'a)|'' and Aa = [-61a, tt - 
Oa] \ H\(Ja\ - aa)/M,TT{aa " Wa\)/M] . We have 

/ h{x)dx= [ H{x)dx, [ \h'{x)\dx= [ \H'{x)\dx. 

J Da JAa JDa J Aa 

and by the 7r-periodicity of the integrand H we have 



\ H{x)dx= / H{x)dx+ / H{x)dx 

JAa J -9a JlTSa/M 

= / H{x)dx+ / H{x)dx = 

Jit— 9 a Jnsa/M Jtt 



TT-TTS^/M 



'TTSa/M 

Analogously, 

/'TT — TTS /M 

/ \H'{x)\ dx= " \H'{x)\dx. 

JAa JTrSa/M 

For X G [TTSa/M, TT — TTSa/M] the sine is positive and 

\H'{x)\< |5- 2|sin9-3(a;)|cos(a;)| -|-gsin«-2(a;) 
It is easy to check that for g ^ 2 we have 

fn—ns^/M 

l-KSa/M 

\q-2\[ I sin«-^(a;)dsin(a;) 



H{x) dx. 



f- 



g-2|sin« ^ {x) \ cos{x) \ dx 

/2 



v'Sa/M 



■K — TTS /M 



Tr/2 



sin* ^(a;) dsin(a;) 
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Prom this we get 



/•TT — TTS /M / / — 

~^ \H'{x)\ dx<{l- 5,,2) (2 + sin''-^ [^-^ 



+ sin' 



g-2 [ 

M 



I [ sin«-2(a;) 
Jo 



dx. 



We then finally get 

f^-K — irSa/M 



M 



J-nsJM 



dx 



<^ (l-^,,2)(2 + sin«-^ 



M 



sm 



q-2 ( ^ 
M 



■ sm 



q-2 ^ 



Observe also that 



sin(7rSa) = sin(7rSa) = sin(7rSa)- 
Since sin(a;)/[Msin(a;/M)] < 1 for a; e (0,7r], we get 



7rsin(7rSa) 



M 



SM,q - sin(7rsa) 



TT — TT S ^ / M 



H{x) dx 



< 7r(l - (5g,2) sin«-i 



M 



■ sm 



9-1 



M 



TT sm^ 



M 



+ 



7rsin(7rSa) 
M 



2(l-(5g,2) + «^ sin«-2(a;)(ia;j. 



Using sin(7rSa/M) < tt/M we obtain 



7rsin(7rSo) 



M 



SM,q - sin(7rsa) 



7T — TTS_^/M 



TT'Sa/M 



H{x) dx 



^{1 - 5g,2) + q sm"-^{x)dx 



Finally, since sm'^{Mda) = sin^(7rSa), we complete the proof by using the 
estimate of SM,q m (1.4). □ 

Theorem 1 implies the following corollary. 

Corollary 1. Let q & (l,oo). IJoa^l then e'^^^if^M) =O.IfaaiZ then 
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^ I J sin«-2 (x) I sin(x + 20a)\ 



+ 



sin(7rsa) 

^min(l,g-l) 



dx 



, (1.6) 



with Sa G (0, and i/ie factor in the big O notation is independent of f 
from Bjv, and also independent of N. 

We now consider the case q = 1 and present estimates of e'^^{f, M) in the 
following lemma. 

Lemma 1. Let a = af. If aa € Z then e^^^{f, M) =O.Ifaa^I' then 



S\T?{lTSa) Sin(20a) 



M 



-'M,a 



< 



sin (TTSg) 
M ' 



cos(2^„)|, (1.7) 



where Sa € (0, |], and 



M-l 



j=0 



cot 



+ Sg) 

M 



Proof. The case (Ja G Z can be proved as in Theorem 1. Assume that Ga ^ 
Using the form of p/(j) from Section 1.2, we have 



2M2 



sin(7r(.7+c7a)/M) 



sin(7r(j - aa)/M) 



sin(7r(,7 - cra)/M) 



sin(7r(j + CTa)/M) 



As in the proof of Theorem 1 we conclude that 



M2 



where 

M-l 

Sm,i = 

i=o 



sin(^(j + aa)/M) 



sin(7r(j - <Ta)/M) 



M-l 

E 

i=o 



sin(7r(j - [aal + Sa)/M + 26'a) 



sin(7r(i- \aa\ +Sa)/M) 



with Sa = [cTa] — (Ja- Changing the index j in the second sum to j — \<Ja\ , and 
using periodicity of the sine, we get 



M-l 



'S'm.i = 

j=0 



sin(7r(j + So)/M + 26'a) 



sin(7r(i + s„)/M) 



and consequently 
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M-l 



Sm,i = 

j=0 



cos(26'a) + sin(26'o) cot 



7r(j + Sg) 

M 



Using the triangle inequality twice, we obtain 

M-l 



Sm,i - sm{20a) ^ 

j=0 



cot 



+ Sg) 

M 



< M|cos(26la) 



Let Sg = (Tg— [ag\ . Observe that Sg = 1 — Sg. Since the cotangent is 7r-periodic 
and the function | cot(7r(-)/M)| is even, we get 



M-l 

E 

j=0 



cot 



7r(j + Sg) 

M 



M-l 

= E 

3=0 



cot 



M 



= MEM,a- 



This and 



sin^(M^a) = sin^(7rcra) = sin^(7rSa) 
yield (1.7) as claimed. 



□ 



Prom Lemma 1 we see that the sum SM,a is the most important part of the 
local average error e*^^(M, /). We now estimate SM,a- 

Lemma 2. Assume that Ug^'L and M > 3. Then 



M 





1 


(^) 


" M 



cot 



7r(M - 1 + Sa) 



M 



_1 p 



TT{M-1 + Sa)/M 



I cotxl dx 



(1 + Sa)/M 



< 



ttM 



TV(M-1+Sa)/M 



■rT(l+Sa)/M 



sin^ X 



dx. (1.8) 



Proof. This can be shown by applying the error formula for rectangle quadra- 
tures (1.5). Note that TrEM,a - jj cot{nSg/M) - fj\ cot(7r(M - 1 + Sg)/M)\ 

is the rectangle quadrature for the integral J^^^^^y^^^'^ \cotx\dx with 
fc = M - 2 > 1. We then obtain (1.8) by using (1.5). " □ 

We now present the final estimate on the local average error e*^^(/, M). 

Theorem 2. Assume that / e Bjv and a = Uf. For M >3, the average error 
of the QS algorithm for the function f satisfies 



2sin^(7rsa)sin(26la) InM 



M 



Stt + 2 + ln(7r^) . , ^ 



Proof For cTq € Z wc have Sg = and (1.9) holds since e^^^if, M) = by [4]. 
Assume that (7^ ^ From Lemmas 1 and 2 we have 
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sin2(7rs„)sin(2e„) /-(M-i+«a)/M 

(\+Sa)/M 



er(/,M)- \cotx\dx 



< 



sin^(7rsa) 
M 



sin(26>a) 
M 



cot 



M 



cot 



7r(M -l + Sa) 



+ 



^ r-K{M-l+Sa)/M ^ 



7i-(i+sa)/M sm X 



M 



■dx + |cos(26'a)| 



Observe that 

fn(M-l + Sa)/M 



L 



I cota;| da; = In I sin 

(l+s„)/M 

7r(M - 1 + s„ 



• -l/7r(l + So)\ . _i/7r(l -Sa) 

— ' ' sm ' 



M 



M 



cot 



M 

(M-l+s„)/M 



7r(l+Sa)/M 



Sin a; 



cot 



dx = cot 



^(1 - Sg) 
M 

7r(l - fig) 

M 



< cot 



+ cot 



7r(l + Sa) 

M 



< 2 cot 

The four formulas above yield 

sin^(7rsa)sin(26'a) 



M 



er(/,M) 



X In sin 



• -1 + Sa)\ . _i/7r(l -Sa) 



Sin 



M 



< 



sin2(7rs„) / (2 + 2/it) sin(26l„) /tts^ 



M 



cot(^^ ) +|cos(2^„) 



Observe that sin(7rSo)/[Af sin(7rSo/M)] < 1 since Sa € (0, 5]. This and the 
obvious estimates of sine and cosine yield 



sin^(7rsa) sin(26'a) 



ttM 



X In sin 



M 



sm 



M 



< 3+- 



2 \ sin(7rsa) 



TT J M 



. (1.10) 



Consider now the left hand side of (1.10). Remembering that M > 3, and 
since 2x/Tr < sina; < x for a; G [0, we get 
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<ln(.^). (1.11) 

Thus by (1.10) and (1.11) we get the final estimate (1.9). □ 



1.3.2 Worst- Average Error 

Prom Corollary 1 and Theorem 2 we get sharp estimates on the worst-average 
error of the QS algorithm. 

Theorem 3. Let M > 3. Then the worst-average error of the QS algorithm 
satisfies the following bounds. 

• For q € (1, oo), 

^(M) < ^(ll\ui^-^{x)dx^ ^\l + o{l)). (1.12) 
The last estimate is sharp, i.e., 

e~^(M)=©(^). (1.13) 
In particular, for M — 2 divisible by 4 we have 
e~*^(M)>^(^i^\in'^-^(x)|cos(x)|''da;) ^\l + o{l)), (1.14) 

and the ratio of the integrals in (1-12) and (1.14) o.^^ approximately 1 for 
q close to 1. 

• For q = 1, 

This estimate is sharp, i.e., 

gwor-avg(^) = 0(M-i InM). (1.16) 
In particular, for M — 2 divisible by 4 we have 

wor-avg , .... 2 In M 37r + 2 + ln(7r^) 
^ > - n M Mtt 

Proof. Consider first the case q G (l,oo). By Corollary 1 we have for all 



gavg 



if,M) < J\u.^-^x)dxy\l + o{l)) 
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where o(l) is independent of /. This yields (1.12). 

The estimate (1.12) is sharp since we can take a Boolean function / such 
that « |. Then (1.6) yields (1.13). In particular, for M = 4A: + 2 and 
Qf = 1/2 we have 9af = 7r/4, = M/A = fc + 1/2 and Sa = |. Therefore 

^^ = W7^{1 [ ^^^''~\^) I cos(a;)r dx^ '\l + 0(1)) 

which proves (1.14). For q close to 1, the value of /J^ sin^^^(x)(ix is mostly due 
to the integrand values close to and tt. Since | cos(a;)|'^ is then approximately 
equal to one, the ratio of the upper and lower bound integrals is about 1. 

For q = 1 the estimate (1-15) follows directly from Theorem 2. To 
prove (1.16) it is enough to choose a Boolean / for which the numbers 

sin2(7rsa^)sin(26'„^) = sin2(Af6»a^) sin(20„^) 

are uniformly (in M) separated from 0, see Theorem 2. More precisely, since 
a/ can take any value k/N for fc = 0, 1, . . . , TV, wc take a Boolean function / 
such that |a/ - sin^(7r/4 + 7r/(5Af))| < l/(2iV). For sufficiently large N, we 
have ~ jtt + ^^tt. For large M = Ak + (3 with (3 e {0, 1, 2, 3}, we then 
have 

sin2(M^„^)sin(2e„^) « sin^ (^^tt^ sin (^1^ + ^^n^ > c> 0, 

for some c independent of M. 

In particular, for M — 2 divisible by 4 we take N > M and a Boolean 
function / € Bjv with a/ = 1/2. Then 

Sa = I and sin^(7rsa)sin(26la) = 1, 

which leads the last estimate of Theorem 3. □ 



1.3.3 Quantum Summation Algorithm with Repetitions 

The success probability of the QS algorithm is increased by repeating it sev- 
eral times and taking the median of the outputs as the final output, see e.g., [2]. 
We show in this section that this procedure also leads to an improvement of 
the worst-average error estimate. 

We perform 2n-t- 1 repetitions of the QS algorithm for some n G {0, 1, . . .}. 
We obtain sin^(7rji/M), sin^(7rj2/M), . . . , sin^(7rj2n-i-i/Af) and let a„ j be the 
median of the obtained outputs, i.e., the {n + l)st number in the ordered 
sequence. Let pnj{a), a € Am, be the probability that the median a„j is 
equal to a. This probability depends on the distribution function Ff of the 
original outputs from Am = {sm^{Trj /M) : j = 0, 1, ... ,M — 1}, which is 
defined as 
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Ff(a) = /^"'e^M,a'<a Pfi^') for a > 0, 
^ \0 fora = 0. 

It is known, see [8] p. 410, that the distribution of the median a„ j is of the 
form 

= (2n + l) ' / Vae^M. (1.17) 

Wc arc now ready to estimate the worst-average error of the QS algorithm 
with 2n + 1 repetitions 



e~HM) = sup max ( ^ p„,/(a)|a/ - aA '\ q € [l,oo). 



We estimate e™''~''^s(M) by using Theorem 12 of [1] which states that 
the QS algorithm with M queries computes a/ such that 

— a/I > ''^TyJ with probability at most — 

for any positive integer k, Here ci and C2 are absolute constants and / is any 
Boolean function from Bjv- If 

k 

\af-anj\ > ci— , 

then for at least n outcomes af{ji), . . . ,af{jn) we must have 

k 

\af -af{ji)\ > ci— for / = l,...,n. 

But the probability that this occurs is bounded by 

2n+ 1\ /C2\" 



(I) ■ 



n 

It follows then that (with c which may depend on n) 

Prob{ |a/ - a„,/| > ak/M} < ck~"-. 
We now use the standard summation by parts. Define 

Pk = Prob{ci(A;- 1)/M < |a/ - a„j| < Cik/M}. 
Then, by the above estimate we get for an arbitrary integer I, 



k>l 
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Therefore 



e--— s(M)« < ^pfc(cifc/M)9 = (ci/M)« ^p, - (/ - 1)") 

oc oo 

1=1 k=l 

oo 

1=1 

for n > g and with the number c = Cq,„ depending only on q and n. In fact, 
taking n = [g] + 1 it is easy to check that Cg_„ is a single exponential function 
of q. Hence, by taking the gth root wc have 

with Cgf',? of order 1. Therefore we have proven the following theorem. 

Theorem 4. The worst-average error of the median of2{\q'] + 1) + 1 repeti- 
tions of the QS algorithm with M quantum queries satisfies 

with an absolute constant in the big O notation independent of q and M. □ 

The essence of Theorem 4 is that the number of repetitions of the QS 
algorithm is independent of M and depends only linearly on q. Still, it allows 
to essentially improve the worst-average error of the QS algorithm. As we 
already mentioned in the introduction, the bound of order is a lower 

bound on the worst-average error of any quantum algorithm. Hence, the QS 
algorithm with repetitions enjoys optimality also in the worst-average setting. 
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